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We show analytically that the heat conductivity of oscillator chains diverges with system size N as 
iV 1 / 3 , which is the same as for one-dimensional fluids. For long cylinders, we use the hydrodynamic 
equations for a crystal in one dimension. This is appropriate for stiff systems such as nanotubes, 
where the eventual crossover to a fluid only sets in at unrealistically large TV. Despite the extra 
equation compared to a fluid, the scaling of the heat conductivity is unchanged. For strictly one- 
dimensional chains, we show that the dynamic equations are those of a fluid at all length scales 
even if the static order extends to very large N. The discrepancy between our results and numerical 
simulations on Fermi-Pasta-Ulam chains is discussed. 

PACS numbers: 



I. INTRODUCTION 

When a small temperature gradient is applied across a 
system, the heat current that flows through it is expected 
to follow Fourier's law, j = —nVT However, for one- 
dimensional systems there are numerous l3l examples of 
systems where numerical simulations H, IE I B , El J5 p. fiol. 
Elm HQ or exact analytical results pUlrl fl7piave 
shown that Fourier's law breaks down. The thermal con- 
ductivity k is observed to be anomalous, i.e. dependent 
on the system size and diver gen t in the thermodynamic 
limit. Despite the existence [T7j of an analytical renor- 
malization group (RG) prediction for the behavior of k, 
there is still considerable controversy about whether the 
theory applies to chains of stiff oscillators, where the par- 
ticles form a lattice to an excellent approximation. Since 
this is the appropriate description for nanotubes, this 
class of systems is of great topical interest. An under- 
standing of k also has broader implications, since other 
transport coefficients can also be anomalous in low di- 
mensional systems [l7|. 

Numerical_|i|l 11 |1J1M El 13 13 13 and ana- 
lytical [3 ElE! 13 EI El] studies show k oc N a for 
various one-dimensional systems, where N is the size of 
the system and a > 0. A renormalization group (RG) 
analysis |17| has shown that for one-dimensional non- 
integrable systems with momentum conservation, a has 
the universal value of 1/3, and that without momentum 
conservation the conductivity is normal (i.e. a = 0). 
This result was derived using the hydrodynamic equa- 
tions of a normal fluid. Simulations on one-dimensional 
hard partic le g ases have yielded a wide range of values 
for a H0 Il6| . This is because very large system sizes 
are needed before the large N behavior is reached, pre- 
sumably because of the limited phase space accessed in 
one-dimensional elastic collisions. This problem has been 
circumvented by using a random collision (RC) model 
where convergence to the asymptotic regime is seen for 
small N. Recent simulations measuring energy diffusion 
in hard particle systems have confirmed that a = 1/3. 

In marked contrast to these results, chains of nonlinear 
oscillators such as Fermi-Pasta-Ulam (FPU) chains [2(j 



and their variants consistently show numerical results 
that disagree with the RG prediction [3 13 13 111 • 
Varying parameters within a given model gives differ- 
ent values for a flCM . but they generally lie between 
0.37 and 0.44 H, li d 111! [3| with occasional exceptions. 
This is contrary to the RG prediction, but is compati- 
ble with the result from a mode-coupling analysis that 
yields a = 2/5 [3 El- 

These numerical results have led 
to suggestions that there are (at least) two universality 
classes for heat conduction in one-dimensional momen- 
tum conserving systems ^2 13 > and that the "fluid- like" 
RG analysis is not appropriate for oscillator chains, which 
are better described as a crystal with small fluctuations 
in particle positions. 

Although any one-dimensional system will eventually 
renormalize to fluid like behavior at sufficiently long 
length scales, for stiff oscillator chains the crystalline 
order persists to very large N. Since the hydrodynamic 
equations for a crystal are different from those for a fluid, 
the large N behavior for such chains could be different 
from the fluid prediction. 

In this paper, we show that this is not the case: a 
should be 1/3 even for stiff oscillator chains. In this 
light, the numerical results on FPU chains are proba- 
bly due to the well-known difficulties with equilibrating 
them [3 HU i necessitating extremely large N before the 
asymptotic behavior is seen. 

In Section II, we first argue that the earlier RG anal- 
ysis E3 applies to chains as well as gases. This result is 
surprising in view of the fact that the standard hydrody- 
namic description of crystals introduces extra degrees of 
freedom for broken symmetries [3 H3 , and is peculiar to 
one-dimensional chains with single-well interparticle in- 
teraction potentials. We then derive the hydrodynamic 
equations for chains with multiple-well interaction po- 
tentials and for quasi one-dimensional systems such as 
tubes [24J. There are four equations instead of three for 
a fluid, because of the extra degree of freedom corre- 
sponding to the broken symmetry. We show that these 
modified hydrodynamic equations also yield a = 1/3. 
In Section III, we elucidate the discrepancy between nu- 
merical results on FPU chains and the analytical predic- 
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tion by performing simulations on FPU chains with col- 
lisions, which interpolate continuously between the RC 
model and pure (collisionless) FPU chains. As we tune 
parameters to move towards pure FPU chains, a new in- 
termediate length scale regime emerges and increases in 
range, while leaving the large N scaling of n unchanged. 
This corroborates the analytical prediction that a = 1/3 
asymptotically, but shows that this will be seen at un- 
usually large N for pure FPU chains. Finally, in Section 
IV, we discuss the applicability of these predictions to 
physical systems. 

II. 1-D HYDRODYNAMICS 

The RG analysis discussed earlier ^j) uses the hydro- 
dynamic equations for a one-dimensional fluid |22j, |25| , 

d t p + d x g = 0, 
d t g + d x (gv) = -d x p + -qdlv, 
d t e + d x [(e + p)v] = KdlT + r][{d x vf +vdlv] (1) 

where p, g and e are the local mass, momentum and en- 
ergy densities respectively. The thermodynamic fields, 
temperature, pressure, and velocity are represented by 
T, p, and v respectively. There are two transport co- 
efficients, the viscosity 77 and the thermal conductivity 

K. 

The RG analysis adds noise terms to the second and 
third of Eqs.Q) (the first is an exact identity and there- 
fore has no noise term). Using the Green-Kubo rela- 
tions 26, 27], the critical exponents are derived from the 
equilibrium fluctuations of a system at constant tempera- 
ture. This can be done without approximations by invok- 
ing Galilean invariance and the condition that equal time 
correlations must obey equilibrium statistical mechanics. 

For a crystal, broken symmetries introduce additional 
long-lived hydrodynamic modes [2^, • For a d di- 

mensional crystal, the extra d broken symmetry degrees 
of freedom are associated with |23j the d — 1 transverse 
sound modes, which become propagating instead of dif- 
fusing, and an extra vacancy diffusion mode 23]. This 
vacancy diffusion mode is the only one that survives for 
d = 1. The number of hydrodynamic equations is in- 
creased from three to four. We are thus led to enquire 
whether the results of Ref. 01 are altered for a stiff chain, 
which can be viewed as a one-dimensional crystal. 

Before doing this, we show that for one-dimensional 
monoatomic systems, if the interparticle interaction has 
only one minimum, vacancies do not exist. The fluid 
description of Ref. |17| is therefore valid even if the inter- 
actions are sufficiently stiff that the particles are almost 
on a perfect lattice. This is true for FPU chains. 

This result can be shown on physical grounds: if we 
remove a single particle from an otherwise perfect d di- 
mensional crystal, there is a slight distortion in its vicin- 
ity, but a vacancy remains. This vacancy only moves 
when a neighboring particle hops into the empty site, re- 
sulting in very slow diffusive motion of vacancies (and 



interstitials) through the system. In contrast, if a parti- 
cle is removed from a FPU-like one-dimensional chain, its 
neighbors move immediately, on microscopic — i.e. non- 
hydrodynamic — time scales to close the gap. Although 
the reduction in mass density caused by removing the 
particle is not eliminated, but only smeared out, this 
is what one would expect, since mass density is con- 
served in Eqs.Q. The extra degree of freedom intro- 
duced in the standard hydrodynamic treatment of crys- 
tals [52, H3, HE , distortions of the broken symmetry, 
is thus indistinguishable from density fluctuations, and 
the fluid description of Eqs.Q) remains valid. 

There are two ways in which the argument of the previ- 
ous paragraph can break down. First, if the interparticle 
potential has multiple minima, it is possible to create a 
situation where the separation between adjacent parti- 
cles is Xi + \ — Xi — a + 5i_ n b for some n, and the force 
on all the particles is zero. Here a is the lattice constant 
and b is the separation between minima of the interpar- 
ticle potential, in the simplest case equal to a. The re- 
sulting vacancy at site n will diffuse slowly through the 
system and is therefore a hydrodynamic mode, distinct 
from lattice distortions and momentum which combine 
to produce longitudinal sound. Second, if the system is 
only quasi one-dimensional, for example a cylinder, va- 
cancies can be introduced as for d > 1 . Neither of these 
exceptions applies to FPU chains. 

The second exception is, however, applicable to nan- 
otubes. In view of their importance, we derive the hydro- 
dynamic equations for a one-dimensional crystal, modi- 
fying Eqs.QJ, and show that the heat conductivity expo- 
nent a is still not changed. This is because the assump- 
tions that allowed a to be calculated exactly for a fluid 
are still valid for the crystal. (Nanotubes have additional 
transverse and torsional degrees of freedom, but we will 
argue later that these do not affect a either.) 

A hydrodynamic theory is formed from the continuity 
equation of the conserved quantities and the time deriva- 
tive of all the broken symmetries. In the case of the 1-d 
crystal, the equations are 

dtp + d x g = 0, 
d t g + d x n = 0, 
d t e + d x f = 0, 
d t u = v+j u , (2) 

where u is the displacment and the currents, 7r, j e , and 
j u , are still to be determined. Galilean invariance de- 
mands the inclusion of the v in the last equation. 

Following the methods from j2f| and |2g , constitutive 
equations can be derived from entropy arguments. The 
first law of thermodynamics for this sytem can be rewrit- 
ten in terms of densities, 

Tds — de — fidp — vdg — hd(d x u), (3) 

where s is the entropy density and p is the chemical po- 
tential per unit mass. A uniform translation cannot alter 
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the energy, therefore energy can only depend on gradi- 
ents of u. The conjugate field to d x u is h = dE/d(d x u) 
at constant s, p, and g. The displacement variable also 
changes the pressure to 

p = —e + Ts + up + vg + hd x u. (4) 

Using these two relations and the continuity equations, 
we can derive an entropy "continuity" equation with 
source terms, 

T[d t s + d x (vs + 9.)] = ~d x T -(g- pv)d x p 

+ (j u + vd x u)d x h - (n + h-p- gv)d x v, (5) 
Q = j e — p(g — pv) — irv + hj u — ev + gv 2 + hvd x u. (6) 

By integrating both sides of Eq.JSJ over a large volume 
where the entropy "current" vs + Q/T is negligible at the 
surface, we have an expression for the time derivative of 
the total entropy, 

+(j u + vd x u)d x h -(n + h-p- gv)d x v}. (7) 

The constitutive equations for the currents can now 
be derived by using the condition of entropy creation, 
dS/dt > |23| . The equality holds only when there is 
no dissipation. Without dissipation, each term of Eq.JTJ) 
must independently be zero, giving the reactive terms of 
the currents: 

.9 = pv, 
j% = -vd x u, 
7Tr = p-h + gv, 
Qr = 0, 

j% = (e + p-h)v. (8) 

The —vd x u term in the displacement current does not ap- 
pear in the derivation of crystal hydrodynamics from |25j 
but appears to us to be correct. With dissipation, en- 
tropy is produced and the currents have dissipative com- 
ponents: 

f D = Td x h-Xd x T, 

7i\D = -rjdxV, 

Q D = -nd x T-XTd x h, 

f D = -(K-Xh)d x T-(XT + Th)d x h-T]vd x v. (9) 

The transport coefficients k, r), A, and T are all positive. 
In general, there are higher order terms in the gradient 
expansion, but we will only retain these first order terms 
(which is valid for small gradients). The mass current g, 
does not have a dissipative term because it is a conserved 
density itself. 

The displacement field introduces two new transport 
coefficents: one associated to the relaxation of its conju- 
gate field h and the second due to the possibility of cross 



coupling between the heat current and the displacement 
field (and the displacement current and temperature 
field). This cross coupling is possible because the field 
and the current have opposite signs under time-reversal, 
a necessary condition for dissipative currents 22] . For 
this reason, there are no cross couplings with n. Using 
these expressions for the currents in the continuity equa- 
tions closes our set of equations, 

dtp + d x g = 0, 
d t g + d x {gv) = -d x p + d x h + r]d x v + ( g , 
d t e + d x [(e+p- h)v] = nd 2 x T + V [{d x v) 2 + vd 2 x v] 

+v[{d x h) 2 + hd 2 x h] + x[Td 2 x h - hd 2 x T] + Ce, 

d t u + (d x u - l)v = Td x h - Xd x T + Cu- (10) 

Eas- HlUfl are the full nonlinear hydrodynamic equations 
of the 1-d crystal. Noise, represented by the £'s, is intro- 
duced to the last three equations for the RG analysis. 

The hydrodynamic modes and their dispersion rela- 
tions can be solved for the linearized theory. The three 
variable normal fluid has three modes, a heat diffusion 
mode and two propagating sound modes . The addi- 
tional displacement variable now introduces another dif- 
fusive mode representing vacancy diffusion. A vacancy 
diffusion mode is seen in the longitudinal hydrodynamics 
of the three dimensional crystal as well [23L I2M |28[ . In 
these linearized theories, quantities such as the suscep- 
tibilities and correlations can be calculated and Green- 
Kubo relations can be derived • Although the details 
of such calculations are different with four hydrodynamic 
modes instead of three, Eas. (|10|) are clearly still Galilean 
invariant, and — when no temperature gradient is im- 
posed externally — have equal time fluctuations that are 
drawn from the canonical ensemble. As in Ref. the 
various critical exponents are therefore determined and, 
in particular, a = 1/3. 

In summary, for one-dimensional momentum- 
conserving crystals (that reach local thermal equi- 
librium), Eqs.Q describe the hydrodynamics when 
vacancies are not possible, as is the case for FPU chains. 
When vacancies are possible, Eas. (|10|) apply, but the 
symmetry arguments that determine a = 1/3 for Eqs.JQ) 
are still valid. In the next section, we numerically show 
the possible source of the discrepancy between prior 
numerical results on FPU-like lattices 0, 0, 0, and 
this exact result. 



III. NON-EQUILIBRIUM SIMULATIONS 

Results of past numerical simulations [il fiol ITU [T3. 

consistently measure an exponent larger than 1/3 for 
FPU- like chains. Due to the exact derivation of the 
analytical prediction, the good agreement between this 
prediction and gas simulations 0, B 1 an d the well- 
known problems of convergence in one-dimensional sys- 
tems [23, we conjecture that the FPU simulations have 
not reached the asymptotic large N regime. Going to 
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larger and larger system sizes to confirm this would 
be prohibitive. Instead, we use a tunable model which 
smoothly interpolates between the RC gas and the FPU 
chain, where this is easier to see. 

We introduce FPU-like springs to the random collision 
model. As the springs become stiffer, the system contin- 
uously evolves from the pure RC gas with a — 1/3 to the 
collisionless FPU chain: from a fluid to a "crystal". As 
the spring constants are increased, we observe that an 
intermediate length scale regime emerges, whereas the 
large N scaling of K remains ~ A 1 / 3 . 

For the simulations, we numerically integrate a system 
of N particles with nearest neighbor interactions. The 
interaction potential between every pair of particles is of 
the generalized FPU form 

n*) = ^+ 1* 8 (ii) 

where z — Xj+i — Xi is the compression (or elongation) of 
the spring and Xi is the deviation from the equilibrium 
position of the i th particle. Collisions can be neglected 
in the limit when the lattice constant a is large, but for 
finite a and finite temperature, the probability of two 
particles colliding is nonzero. When two particles come 
into contact, z = — a, they undergo a random collision as 
in Ref. 0. 

For the simulations, we only use terms up to fourth or- 
der in Eq. i|ll|) . In their original work, Fermi, Pasta, and 
Ulam 20] noticed that the odd frequency normal modes 
do not mix with the even frequency normal modes when 
the interparticle potential is even. To avoid this problem 
and any others that may arise from this accidental sym- 
metry of an FPU-/3 chain, we use a finite k% in all of our 
chain simulations. 

Heat baths are connected to the end particles to main- 
tain a temperature gradient across the system. The scal- 
ing of the conductivity is determined by measuring the 
size dependence of the current maintaining a (small) tem- 
perature difference between the baths. For our numerics, 
the leftmost (i = 1) and rightmost (i = N) particles are 
connected to Nose- Hoover heat baths [3(3, [jy} . The aux- 
iliary degree of freedom introduced in the RC model, the 
transverse momentum, does not couple to the baths or 
enter the dynamics except in collisions, and so is not seen 
in the pure FPU (collisionless) limit. 

We adopt a convention in which the energy contained 
in each spring is symmetrically shared by the two par- 
ticles attached to it. The first and last springs are only 
attached to one particle and a wall, therefore their en- 
tire energy is attributed to the attached particle. Using 
this convention for the energy density, the energy cur- 
rent consists of two parts. The first is an advective part, 
j adv (x,t) = Y^,i eiViS(x — Xi(tj). The second is piecewise 
constant and jumps at each particle with ji+n — ji,i-i — 
-e l . From this, j i+lti = -(l/2)(v i+ i + Vi)V'(xi+i - afj), 
with jo,! and jN,N+i equal to the energy currents from 
the reservoirs to the end particles. Thus the total current 



flowing through the system is 

JV - N 

jN = ^e l v i --^(x l+1 -x l +a){vi + i+v i )V'(x i+ i-x i ). 

i=l i=l 

(12) 

It is possible to show that the time average < jN > is 
equal to -(a/2) ^(Vi+i + Vi)V'(x i+ \ - x t ). However, 
with collisions it is more useful to keep Eq. (|12fl in its en- 
tirety, since 2^+1 — Xi+a = at collisions, where V is sin- 
gular. In fact, for hard sphere particles, jN = e i v i- 

Numerical integration of this system is complicated 
by the combination of collisions and springs. Without 
springs, the particles travel freely and a fast event driven 
simulation can be implemented For a collisionless 
system, the equations of motion can simply be integrated 
numerically by standard algorithms. For our simulations, 
the fourth order Runge-Kutta algorithm is used. Because 
of the occurrence of collisions, our algorithm checks for 
collisions after every trial timestep. If a collision is seen 
to have occurred, the trial step is discarded, and the in- 
tegrator evolves the system by a smaller stepsize to the 
(extrapolated) point at which the two particles collide. 
Because of the need for trial steps, integrating systems 
with collisions and springs requires substantial amounts 
of computational time. 

We use a stepsize of dt = 0.01 for all simulations and 
check that the results do not change for smaller stepsizes. 
We allow the system to equilibrate for ~ 10 s steps before 
measurements are taken. Each measurement consists of 
the particle and time averaged current, where a block 
of 10 s timesteps is used for the time average and 10 5 
steps separate each block. Each data point in figure |3 
is the mean of many such measurements and the error 
bars shown are the root mean square errors of the block 
measurements. 

The parameters in the model are the masses of the par- 
ticles nii , the lattice constant a, the temperatures of the 
baths Ti, and Tr, and the spring constants k n in Eci. (|llf) . 
Alternating masses are used due to the fast convergence 
seen in the pure RC model for such systems ||. In par- 
ticular, we use a mass ratio of 2.62 and nij = 1 for all 
odd i. The bath temperatures are 1.2 and 1.0 for the 
left and right baths, respectively. For these tempera- 
tures, we have checked that the system is in the linear 
response regime. With the temperatures of the baths 
and the masses of the particles fixed at these values, the 
collision rate can be altered by changing the lattice con- 
stant and/or spring constants. We have chosen to fix the 
lattice constant to a — 1/2. 

Figure ^ shows the kinetic temperature profile after 
steady state is reached for systems of size N — 256, 
showing increasing curvature as collisions decrease. For 
the pure FPU system, the masses of all the particles are 
chosen to be the same. This is because with alternating 
masses, the kinetic temperature for particles in the pure 
FPU chain also shows an odd-even oscillation, indicating 
that equilibrating the system with Nose-Hoover baths is 
problematic. All other systems simulated have alternat- 
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FIG. 1: The temperature profile for a pure RC gas, a mixed 
system with springs and collisions, and a collisionless chain 
with springs. The system size N is 256. The spring constants 
(fc 2 , fcs, ki) of Eq.JTTJ are (I, 0.1 , I). 
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FIG. 2: Log-log plot of (jN) versus N, where j is the time- 
averaged heat current and N is the number of particles, for 
systems with different spring strengths. The spring constants 
(k2, kz, ki) of Eq. IIH are specified in the key. The figure also 
shows a collisionless FPU chain at the top for comparison. 
The systems without springs or weak springs have slopes of 
~ 1/3 at large N. The slope increases with spring strength 
and is > 0.4 for the pure FPU system at the largest N shown. 
When shown, the error bars are root mean square errors. The 
plots are shifted for clarity. 



ing masses; as mentioned in the previous paragraph, this 
allows the large N limit to be reached faster. We verify- 
that the temperature profiles for these do not oscillate, 
as shown in Figure ^ 

After the temperature gradient is established, the en- 
ergy current is measured. Figure [21 shows a log-log plot 
of the current times the number of particles (jN) versus 
N; the slope of the plot determines a. 

For no springs, the pure RC model has a slope of ~ 1/3, 
similar to the earlier results in 9] . This asymptotic value 



is reached for relatively small systems (N > 128). As 
the springs are turned on, an intermediate length scale 
regime emerges and the asymptotic large N regime moves 
out, but the slope in the large N regime is still 1/3. For 
the strongest springs shown in Figure [2 the asymptotic 
regime has moved out of the range simulated. Figure [21 
also shows the results for a pure FPU chain, which agrees 
with earlier results [H El |H 13. The slope of > 0.4 
for this system is comparable to the strong spring system 
with collisions, supporting the assertion that this is an 
intermediate length scale phenomenon. As discussed in 
the previous paragraph, the masses of all the particles in 
the pure FPU chain are chosen to be equal. This is the 
cause for the different behavior for very small N, where 
the slope increases with N instead of decreasing. 



IV. DISCUSSION 

In the light of the experience with hard particle gases 
and the random collision model, it is natural to expect 
that introducing extra degrees of freedom would cause 
FPU chains to show their large N behavior more rapidly. 
In this context, Wang and Li 0] have shown that when 
one transverse degree of freedom is added to the FPU 
chain, the conductivity exponent is numerically seen to 
be 1/3. It would be interesting to see whether other ex- 
tensions of the FPU model can achieve the same result. 
More ambitiously, it would be very useful if a criterion 
could be obtained that would easily determine numeri- 
cally whether a system reaches local thermal equilibrium, 
allowing one to decide whether the hydrodynamic theory 
is applicable. 

The study of low dimensional transport is not merely 
of theoretical interest; the vast potential applications of 
real quasi one-dimensional systems, namely nanotubes, 
demand an understanding of their physical properties. 
Specifically, research and production of electronic devices 
at the micrometer and nanometer scale have shown that 
carbon nanotubes may be used in very efficient cooling 
systems |32j . Experimentally, the thermal conducti vity 
of nanotubes has been found to be extremely high [33] 
which makes nanotubes obvious candidates for cooling 
devices in microelectronics. This potential for applica- 
tions motivates a theoretical investigation of the heat 
transport properties of low dimensional systems. Exper- 
imental [33, 13 an d numerical [3f|, |3(| studies of nan- 
otubes have shown that the phonon contribution to heat 
transport is much larger than the electronic contribution, 
so that the fact that we have only considered lattice mo- 
tion in this paper is not a problem. 

The two dimensional nature of the tube allows for 
the existence of vacancies. The tubular shape also in- 
troduces other transverse hydrodynamic modes (a bro- 
ken symmetry and conserved momentum for torsional 
and the two transverse motions), but as mentioned ear- 
lier, these modes do not affect the symmetry arguments, 
and the heat conductivity exponent is expected to still 
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be 1/3. The longitudinal hydrodynamics for this quasi 
one-dimensional system is the four component theory, 
Eqs.Jini. 

There are a few concerns in applying the prediction 
of a — 1/3. First, on short length scales, the phonon 
motion is ballistic. The phonon mean free path in car- 
bon nanotubes is ~ 1/im |34j . and our discussion in this 
paper would apply to systems that are larger than this 
scale. Recent thermal conductivity measurements have 
only been made on systems with lengths in this order of 
magnitude |34l | . The unusually high conductivity in such 
systems is correctly attributed to the ballistic transport 
nature at these scales [34]] • There have yet to be exten- 
sive experiments on the length dependence of k to our 
knowledge. To test our prediction and measure the ac- 
tual heat conductivity exponent, length dependent exper- 
iments of long (3> Ifim) nanotubes must be done. Sec- 
ond, even beyond this length scale it is possible (based 
on what we have seen for FPU chains) that equilibration 
is imperfect. However, this is less likely to be a problem 
with the extra transverse modes [bH ] which are present 
in a nanotube, and in any case, would only change a 
slightly. Third, under renormalization, the system flows 
to its fixed point only for large N. It is not clear what 
limitation this imposes, but the nonlinearity in the hy- 
drodynamic equations that causes anomalous scaling be- 
havior is the advective term, whose strength is unity and 



whose effect should therefore be seen even for small sys- 
tem sizes j^E |H3 ■ 

In this paper, we have shown analytically that one- 
dimensional chains of particles connected with nonlinear 
springs have a heat conductivity that diverges as a func- 
tion of the chain length N as ~ TV 1 / 3 , which is the same as 
the earlier result for hard sphere particles and other fluid 
systems. For quasi one-dimensional systems, this result is 
obtained from the crystalline hydrodynamic equations. It 
does not rely on the eventual crossover from a crystal to a 
fluid which must happen for any one-dimensional system. 
For Fermi-Pasta-Ulam chains and other examples where 
the intcrparticle potential has a single minimum, we have 
obtained the stronger result that the fluid hydrodynamic 
theory 01 is applicable on all length scales for the dy- 
namics, even when there is excellent static ordering. In 
the light of these analytical results, the numerical results 
on FPU chains that show a heat conductivity exponent 
of ~ 0.4 is probably due to imperfect equilibration. 
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